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The possibility of enhancing the ratio of output to 
input power Q in a simple mirror machine by polariz- 
ing deuterium-tritium (D-T) nuclei is evaluated. Tak- 
ing the Livermore mirror reference design mirror ratio 
of 6.54 , the expected sin 2 d angular distribution of fu- 
sion decay products reduces immediate losses of alpha 
particles to the loss cone by 7.6% and alpha-ion scat- 
tering losses by -50%. Based on these findings , alpha- 
particle confinement times for a polarized plasma should 
therefore be 1.11 times greater than for isotropic nu- 
clei. Coupling this enhanced alpha-particle heating with 
the expected >50% D-T reaction cross section , a cor- 
responding power ratio for polarized nuclei , Q po iahzed> 
is found to be 1.63 times greater than the classical un- 
polarized value Qdassicat • The effects of this increase in 
Q are assessed for the simple mirror. 


INTRODUCTION 

The availability of multiampere beams of nuclear 
polarized atoms 1-3 and their resistance to depolariza- 
tion under fusion reactor conditions 4 " 9 has sparked in- 
terest in the use of fuel polarization to enhance fusion 
reactor performance. For a magnetic mirror device, po- 
larization of deuterium (D) and tritium (T) nuclei par- 
allel to the B field offers two advantages: 

1 . a 50% increase in the effective nuclear cross sec- 
tion 

2. fusion decay products that are emitted predom- 
inantly perpendicular to the magnetic field in a 
sin 2 6 angular distribution. 

Because of these two effects, polarization of injected 
nuclei is expected to increase the ratio of output to input 
power for the mirror. The higher cross section prom- 


ises 50% higher reactivity and 50% more nuclear power 
for the same conditions (or the same power at lower op- 
erating pressures if beta is seriously constrained). Like- 
wise, the sin 2 d angular distribution of generated alpha 
particles should reduce alpha-particle end losses from 
the mirror. Consequently, more of the alpha-particle 
3520 keV should be deposited in the reactant plasma 
leading to a net reduction in the required external heat- 
ing power. 

Enhanced alpha-particle trapping results from 

1 . the reduced number of alpha particles immedi- 
ately “born” into the loss region 

2. the reduced alpha-ion scattering losses of the re- 
maining alpha particles. 

For these alpha-particle loss mechanisms, analytical 
and numerical enhancement factors are derived, and 
the ratio of polarized to isotropic alpha-particle con- 
finement time is found. These results are incorporated 
into an equilibrium fusion energy balance to translate 
the additional alpha-particle heating into a higher ratio 
of output to input power (henceforth referred to as Q). 


ALPHA PARTICLE HEATING 

As larger experimental fusion machines approach 
reactor conditions, alpha-particle heating is an increas- 
ing issue of concern. 10 Alpha particles generated pri- 
marily by the reaction D + T -► 4 He + n carry fully 
one-fifth of the fusion reaction energy. Since this ki- 
netic energy is “attached” to a charged species, some 
portion is magnetically confined to be partitioned 
among the fusion reactants. Alpha-particle studies prin- 
cipally seek to determine the dynamics of energy re- 
laxation within the background plasma and the total 
amount of energy deposited with each species (see 
Ref. 11 for summary). 

The principal alpha-particle energy loss terms are 
assumed here to be 
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1 . the immediate loss of alpha particles born into 
the loss cones 

2. slowing-down collisions with electrons and ions 

3. pitch-angle scattering by collisions into the loss 
cones. 

Anomalous energy losses such as Alfven wave coupling 
will be ignored. Each loss term is considered in turn 
with the ultimate objective being to develop a simula- 
tion model for the energy coupling of the alpha par- 
ticles to the plasma for both an isotropic and sin 2 d 
alpha-particle distribution. 


LOSS CONE ALPHA PARTICLES 


Ignoring ambipolar potentials, the probability of 
immediate alpha-particle loss for an isotropic distribu- 
tion is P {oss = 1 - cos d 0 , where cos 0 0 - { \ - R~ l ) w2 y 
and R is the mirror ratio. Polarized D-T nuclei, on the 
other hand, have an alpha-particle distribution that is 
proportional to sin 2 d 0 and can be written as 



valid to the first order in do/dti y the fraction of unpo- 
larized nuclei, or as a normalized distribution: 


P(Q) = 


— sin 2 d . 

87T 


(1) 


The probability that an alpha particle is born into the 
loss cone for polarized D-T nuclei is, therefore, 

/> /ass =J- sin 3 d 0 d4> dd 0 = I PWdQ, 

o7T Jo Jo Jo 

where dQ is the differential solid angle dU — sin d dd d<t> 
or 

Ptoss = 1 - jcos t?o(sin 2 do + 2) . (2) 

Table I summarizes the polarized-to-isotropic ratio of 
trapped alpha particles: 

polarized trapped fraction 1 , . 2 « , 
unpolarized trapped fraction 2 

which indicates the alpha-particle trapping enhance- 
ment by D-T polarization. As expected, the trapping 
advantage of the distribution increases with loss cone 
size. 


SLOWING DOWN AND PITCH ANGLE SCATTERING 

The dynamics of the slowing-down and scattering 
process is governed by the Fokker-Planck equation, 
which describes each particle as a material point in 
seven-dimensional phase and time space. Following 
the analysis of Sivukhin 12 and Anderson et al., 13 the 


TABLE I 


Polarized-to-Isotropic Ratio of Trapped Alpha Particles 


Mirror 

Ratio 

0 

Isotropic 

Trapped 

Fraction 

Polarized 

Trapped 

Fraction 

Ratio of 
Polarized 
to Isotropic 
Trapping 

14.9 

15 

0.97 

0.999 

1.03 

4.0 

30 

0.87 

0.983 

1.13 

2.0 

45 

0.71 

0.888 

1.25 

1.3 

60 

0.50 

0.690 

1.38 

1.1 

75 

0.26 

0.382 

1.47 


Fokker-Planck equation is written here as an ordinary 
diffusion equation describing the pitch-angle scattering 
and solved numerically by standard forward-difference 
techniques. Anomalous effects of an energy-dependent 
diffusion coefficient will not be considered. 

The following assumptions are made: 

1 . There is a three-component plasma of Maxwell- 
ian electrons and ions with an alpha-particle distribu- 
tion / that is isotropic in physical space but arbitrary 
in velocity space. 

2. The alpha-electron scattering can be neglected 
relative to ion-alpha scattering (a ~ v~ A , which is small 
for high-velocity electrons). 

3. There is a square-well uniform B field that runs 
parallel to the cylindrical mirror axis and rises discon- 
tinuously at the ends (thus allowing the mirror field to 
appear only in the boundary conditions). 

4. Only the high-energy end of the alpha-particle 
distribution (where df/dt = 0 is a valid assumption) is 
considered. 

Using these assumptions, the Fokker-Planck equation 
takes the form: 


1 3 2 1 d 

v 2 dv V 1 + u sin 0 dd 


(vl ± sin d) = S , (3) 


where I, and I ± are the parallel and perpendicular 
components of the flux density in velocity space, re- 
spectively, and S is a source function, i.e., the number 
density of high-energy particles created per second by 
fusion reactions. Neglecting parallel diffusion relative 
to the perpendicular flux components (valid for t>, < 
v a < v e ), Eq. (3) can be written as 


1 

t s v 2 dv 


[(U 3 + u 3 )/] 


K Vc d 
T s v 3 dx 


<!-*»)? 

dx 


= 5 , 


(4) 
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where 


~ (2ir) l/2 — 
3 m { 


1/2 


e n P 


(kT e ) 3/2 


In A 


tides and the corresponding particle and energy losses 
can be evaluated for isotropic and sin 2 1 ? alpha-particle 
distributions. The unnormalized energy spectrum of 
trapped alpha particles is given by 


v 3 = (9 x) l/2 — [Z]v 3 

m a 

t s = slowing-down time 

v c = critical velocity at which electrons and 
ions contribute equally to the slowing- 
down frictional force 

m e , m a = electron and alpha masses, respectively 

n ey T ey e = electron density, temperature, and charge, 
respectively 

In A = Coulomb logarithm. 

Pitch-angle scattering is described by x = v\\/v = cos d y 
where i? is the angle between v and the B axis. The pa- 
rameters k = Z eff /(2A a [Z]) y Z tff = LZfn i /n e , [Z] = 
L(Zfni/Ai)/n e> and A { = mass of the ions (Aj) and 
alpha particles (A a ) complete the left-hand side of 
Eq. (4). The source function, S = S 0 6(v — v a )H(x) 
[where S 0 is a constant related to the fusion rate and 
H(x) will be given explicitly later] closes the Fokker- 
Planck description. 

Integrating Eq. (4) from v~ to v + reduces the delta 
source to the boundary condition: 


/( v ay x) 


SpHOch^Ts 

yj + v c 


Taking the other two boundary conditions as /( v y x = 
±* t ) = 0 for the loss cone space, the problem is well 
posed. Defining new variables 'F and t as 

* = (V 3 + V 3 )f 



* (*, O 

V 3 (t) + V 3 


dx , 


( 7 ) 


and is shown for / in arbitrary units in Figs. 1 and 2 
for both isotropic (ISOT) and polarized (POL) dis- 
tributions at various mirror ratios. The peak in both 
distributions around 0.1 contrasts with the standard 
Maxwellian distribution in that slow alpha particles are 
preferentially scattered into the loss cones. This altered 
distribution can be understood to arise directly from 
the 1/f 4 dependence of the Rutherford collision cross 
section and becomes more pronounced for smaller loss 
cones. 

The particle loss fraction P y normalized with re- 
spect to the total number of generated particles (includ- 
ing those born into the loss cones), is given by 


P(t) = 



['F (jc, 0) — 'F (x‘, t )] dx 



>F (x,0) dx 


Particle loss fraction as a function of energy during 
slowdown for both isotropic and sin 2 d initial distribu- 
tions in a 10-keV D-T plasma is shown in Fig. 3 for mir- 
ror ratios of 2 and 10. 

The energy fraction Q was found in terms of this 
particle loss fraction by 


Q[PU)] 



for 0 < P < 1 


and 


_ K f 1 + Vc/V* 

“ 3 \1 + ^ 


Eq. (4) can be written as 

d* _ d_ 
dt dx 
with boundary condition 


, d ¥ 

(1 - X 2 ) — 

dx 


and 


( 5 ) 


( 6 ) 


♦ (*, 0 ) = S q vIt s H(x) , |*| < x L 


1 * f (±* z .,/) = 0 , t > 0 , 

where H(x) = 1 for an isotropic alpha-particle distri- 
bution and H(x) = 1 - x 2 for a sin 2 d distribution. 

Numerical solution of Eq. (6) by forward-difference 
methods [step-size ratio At/(Ax) 2 < 0.5 for stability] 
gives 'F as a function of x and t. Using this steady dis- 
tribution function, the energy spectrum of trapped par- 


and is shown in Fig. 4 for mirror ratios of 2 and 10. 

A total energy-scattering rate for the isotropic and 
sin 2 d alpha-particle distribution was found to be -6 
and 2%, respectively, for a mirror ratio of 10, and 12 
and 5%, respectively, for a mirror ratio of 2. The re- 
sults for the isotropic alpha-particle case are in good 
agreement with the findings of Anderson et al., 13 re- 
ported as 6 and 10% for R = 10 and R = 2, respectively. 
Based on this analysis, injection of polarized D-T nu- 
clei has the potential to reduce alpha-particle scatter- 
ing losses by 50% or more. However, given that losses 
are responsible for negligible alpha-particle energy leak- 
age, little improvement in overall Q can be expected 
from this reduction. 

ALPHA PARTICLE CONFINEMENT TIMES 

To find an analytic expression for the ratio of po- 
larized to isotropic alpha-particle confinement times, 
a quasi-Maxwellian distribution function is assumed, 
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Fig. I . Effect of end loss on isotropic steady distribution function. 



Alpha Particle Energy, E/E u 

Fig. 2. Effect of end loss on polarized steady distribution function. 


and Eq. (3) is solved to give the alpha-particle loss rates. 
While on first inspection, the assumption of Maxwell- 
ian alpha-particle behavior may seem artificial, Gal- 
braith and Kammash 14 found a < 1 % difference in the 
key parameter for this investigation, alpha-ion frac- 
tional heating, for a Maxwellian over a mirror alpha- 
particle distribution. Therefore, given the complications 


of a full Fokker-Planck calculation, such a simplify- 
ing assumption seems justified for the purposes of this 
investigation. 

Equation (3) translated to space is 
d d 

P — (sin M a ) + sin d — ( p 2 I p ) = qp 2 sin d , (8) 
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(x 0.1) 



Fig. 3. Particle loss fraction versus alpha-particle energy. 


where q is the source density of alpha particles and l p 
and I & are the momentum and angular flux compo- 
nents, respectively. Following Sivukhin, 12 


1 dF dF 

Id - - &dd TT — A?P T“~ + A#F , 
p dd dp 

^ dF 1 dF 

Ip — — Dnn — — Dpd r ApF , 

pp dP p dd p 


and 


D#p — Dp# . 

As before, the distribution function vanishes at the loss 
cone boundary, F = 0 for d = d 0t and by symmetry, 
dF/dd = 0 for d = tt/2. 

Since an exact solution requires knowledge of both 
the diffusion tensor D a p and the coefficient of dynam- 
ical friction A to find the unknown distribution func- 
tion, the diffusion approximation is employed. By 
replacing D a0 and A by known functions of p and d % 
as if Fwere isotropic, D Pd and A 0 vanish. Further- 
more, Dp P = Z>u, D#d — Al* and A P = A and are all 
functions of p only. Equation (8) can then be written as 


/ . „ dF\ 

■ „ d 

, / „ dF \ 

( sin d — - ) 

+ sin d — 

p^\Di AF ) 

V dd) 

dp 

L V dp )\ 


= - qp 2 sin d . (9) 

Since a quasi-Maxwellian distribution function is 
sought. 


F(p,d) = 6(tf)exp^^L j . (10) 

Equation (9) simply reduces to 



which, when coupled with the boundary conditions, 
gives 


F = 


qp 2 sin d 

In . 

D ± sin d 0 


Comparing this result with Eq. (10) gives 


0 = Cln 


sin d 
sin d 0 


or 


F= C In 



where the normalization constant C is found from the 
particle conservation requirement, 


n — 


/ 


F(d>p) dp . 


( 12 ) 


Since an equilibrium alpha-particle loss term is de- 
sired, the number of alpha particles lost should equal 
the number produced, or 
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Fig. 4a. Energy loss fraction versus alpha-particle energy. 



Fig. 4b. Particle and energy loss fractions versus alpha-particle energy. 


C /* T *^0 ci n i9 

N - I qdp , (13) n = 2%C I In — -sintfrfi? 

*^#0 sin $o 

where the vector dp = 2tt/? 2 sin 2 Odd dp. x f ^ 2 exp { ~P \ ^ 

For an isotropic alpha-particle distribution, 9(i?) = 1, Jo \2mT ) 
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or 


C — 


(2irmT) 


3/2 


In ||cot 




— cos i9 0 


both isotropic and polarized alpha-particle distributions 
[Eq. (14)], and the difference in confinement times 
T = n/N is found in the normalization [Eq. (12)]. In- 
stead of 


such that the loss rate (13) is 

,/2 ln A e 4 n 2 3 ln(V2 + 1) - >/2 


N 




S' 

Jdo 


1 . 0 , sin d 

sin d In — — — dd 
sm uq 


-3/2 


,n ( 


^0 1 q 

cos — ) - cos d 0 


cos d 0 


for the isotropic distribution, the polarized integral is 
1 


For a polarized alpha-particle source term, on the 
other hand, q = q Q (p)^{d) y where *(t?) is normal- 
ized to give the same value for q as for the case of iso- 
tropic alpha particles [^(^) = 1 ]: 

* -jt/2 

^ (d)s\n d dd = cos d 0 , 

'*o 


1 - (^)cos 2 1 ?o 


S’ 

Jdr 


*f T 

J o 0 L- 
X sin d dd 


. sin 0 1 , 2 o 2 q \ 

In (cos d — cos 0 O ) 

sin 0 O 6 


Based on this analysis, the ratio of alpha-particle con- 
finement is 


1 


Tpol 

T iso 


1 — ( I )cos 2 d { 


r* x/2 

0 Jo o 


2 , sin 

- In 

3 sin 1 / o 


1 7 o 

(cos 2 d 

6 


cos 2 do 


sin d dd 


n . 0 , sin d 

sm i? In — — — 
0q sm 0 o 


1 - ( 5 )COS 2 d 0 


i: 

]['■(•“ t) 


COS l?o H cos3 00 

6 


In cot 


?)- 


(15) 


COS t?o 


where ♦(«?) = K sin 2 0 gives 


*(tf) = 


sin 2 d 


1 - ( 5 )cos 2 i ? 0 ‘ 
Solving Eq. (11) for F as before gives 


(14) 


F = 


QoP 

D t 


r» cw r T/2 : 

J# 0 sin d' Ja- 


¥(«T)sint V dd* 


and 


F = 


QoP 


D x (l - ( j )cos 2 «? 0 ) 


- In - - (cos 2 d - cos 2 1 ? 0 ) 

3 sm d 0 6 

Comparing this again with the quasi-Maxwellian form 
gives a 0 dependence of 


0 ( 0 ) = 


1 - ( 5 )cos 2 d 0 
sin d 


Figure 5 summarizes these results for various mirror 
ratios. As expected, Eq. (15) approaches 1 in the limit 
as the loss cone goes to \ and 7t/2, Figure 6 uses an 
approximation formula 1 * for alpha-particle slowing- 
down behavior to demonstrate the effect of this en- 
hanced alpha-particle confinement on the deposited 
plasma energy. 

FUSION ENERGY BALANCE 

To assess the impact of this enhanced alpha-particle 
trapping on calculated Q values, an equilibrium ion and 
electron energy balance is performed. Following Kam- 
mash , 16 the ion and electron injection power per unit 
volume is W si and W S€ , respectively. The geometric 
confinement time is r,- and r e for ions and electrons, 
and when they escape, they carry W Li and W u . 

The term W {J is taken here as the energy exchange 
between each species (W^ = 0 for / = j). With these 
definitions, the ion energy balance can be written as 


[3 

x - In 

[2 


I 


— - (cos 2 d - cos 2 d 0 ) 
sin d 0 6 


- (\n,T,) = 0 = W ia + W si + W Li - W ei 
dt 


Because of the previous normalization of the source 
term, the alpha-particle loss rate N is the same for 


where 


W la = \n}{ov)U„ y (keV/m 3 -s) 
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If U ia is that portion of the alpha-particle energy de- 
posited with the ions, then 


Wsi = Si Vi = 



n 2 (ov > 

4 


(keV/m 3 -s) . 


and analogously for including a C 2 to account for re- 
flectivity of synchrotron radiation from the wall with 
a dimensional size parameter 

D = jS 3 / 2 (Lfl ) ,/2 . 


If the source S, = ion loss rate and r s — thermalization 
time, then 



(keV/m 3 *s) . 


Combining these effects, the ion energy balance is 


(av)(2U ia 4 2 Vi - 


_ ( 127 ; - 8 ^) 
3 Ti) A - 


nn 


9.6 x 10" 18 

* p 


T, 



(16) 


Similarly, the electron energy balance can be written 



= ^ + K, -W u - W x - W c , 

where 

n 2 

W ae = — (ov>U ae (keV/m 3 -s) . 

4 


If £/« is that portion of the alpha-particle energy 
(keV/m 3 -s) deposited with the electrons, then 


m. 


\Te 


(av) 


W ei = 1.2 X 10 


-.8 


-3/2 


T 

1 - 4 

Ti 


W u =-n P T e =- TA~ + 


n e n 2 {ov) 


and 


The electron energy balance simplifies to 


{ov) (2U iot + 2 V e — 3T e ) - (127; - 8V e )/r\T e 


4 


9.6 x 10“ 18 

t-3/2 

1 e 



- 4.8 x 10 ~ 21 C 2 Ty /4 (T e + T,) in 



(17) 


Combined, Eqs. (16) and (17) form a nonlinear set with 
roots equal to the equilibrium ion and electron temper- 
ature if the following quantities are specified: 


1. ion and electron injection energies V h V e 


2. burnup fraction f L = (S, - L { )/S — (1 4 [4/ 
<ai>><A7T>]) -1 


3. ratio of ion to electron confinement times (r// 
Te) = 7 

4. bremsstrahlung parameter C\ 

5. synchrotron radiation reflection adjustment C 2 

6. size parameter D . 

Equations (16) and (17) were solved using a two- 
dimensional Newton-Raphson iterative scheme in 
(7;, T e ) space with the cross-section data called from an 
interpolating subroutine for the smoothed Glasstone- 
Lovberg 17 results. For comparison sake, relevant ma- 
chine parameters were taken from the Livermore mirror 
reference design 18 (LMRD) and are summarized in Ta- 
ble II. 

Figures 7 and 8 show equilibrium temperatures ver- 
sus burnup fraction for both polarized and unpolarized 
fuel injection. As expected, polarized nuclei elevate 
both the equilibrium ion and electron temperatures with 
the more substantial energy gain found in electron tem- 
peratures due to preferential alpha-electron energy ex- 
change. The relevant ratio for Q calculations, however, 
is the ratio of ion energies, which is shown versus burnup 
fraction in Fig. 9. 


dT ( 


W x — C\ - n e — - = (3x 10 ~ 2l )C l n 2 T i 
2 at 


1/2 


Q CALCULATIONS 


where C\ is an adjustable parameter to account for 
changing bremsstrahlung, 

tV c = 6 x 10~ 2S C 2 n 2 Ty /A (T e + 7j) 3/2 (l - 

(keV/m 3 -s) , 


A useful parameter for evaluating mirror reactor 
machines is Q , the ratio of fusion power per unit vol- 
ume ( P F ) to the power per unit volume injected and 
trapped in the plasma (P ,„) . In turn, P in is related to 
the beam energy £j„, and the equivalent injected cur- 
rent by the expression P in = E in I in . For steady-state 
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Electron Temperatures 


Ion Temperatures 


• Polarized 
x Isotropic 


Fig. 7. Equilibrium temperature for isotropic and polarized alpha-particle distributions as a function of fractional burnup. 
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Fig. 8. Equilibrium temperature for isotropic and polarized alpha-particle distributions as a function of fractional burnup. 


operation, the input current must equal the rate of par- where n is the ion density and r is the mean confine- 
ticle losses given by ment time, or 


Iou<= f-dV = f-dV, 
J t J m 


P in = E in f-dV=-^~ f n 2 dV . 
J nr (nr) J 
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Fig. 9. Equilibrium ion temperature for isotropic and polarized alpha-particle distributions as a function of fractional burnup. 


Likewise, since 

PF= <^i„ ldV 

for a 50-50 D-T plasma, Q can be written as 


{rn){ov)E F 
4 E in 


Q polarized 
Q isotropic 


( A 7 * ) enhanced^ ) enhanced 

{ HT y classical ( tJU ) classical 


( E enhanced \ enhanced 

E classical ) ( W*) classical 


(18) 


where the right-hand side of Eq. (18) is the multiplica- 
tive factor desired for a particular design. 


neglecting the differences in (ov) and E in for deute- 
rium and tritium. 

Previous mirror feasibility studies 18 have arrived at 
a reference design values for a mirror Q of 1 .20 with 
77% of the gross electrical power recirculated, a num- 
ber substantially below the economically attractive 
range of 3 to 5 for Q . However, given the fact that 
modest improvements in Q can reduce significantly the 
cost of power, various methods for Q enhancement are 
being investigated (see Ref. 19 for summary). In the 
study of enhanced Q devices, the functional dependence 
of the classical Q predicted by the complete Fokker- 
Planck calculations is retained, but Q en hanced is found 
by multiplying Qdassicat by some constant factor. More 
specifically, since numerous investigations have con- 
cluded that for a fixed injection energy and angle, 

<nr> - E 3/2 log i0 R , 

where E is the mean ion energy and R is the effective 
mirror ratio; Qpoi/Qisotropic can be written for cases of 
interest as 


TABLE II 


Reference Mirror Parameters 


Physics Parameters 

Engineering Parameters 

Mirror ratio = 6.54 
P c = 3.52 keV 
P s = 14.1 A/keV 
91 % alpha trapping 
<<n>> a = 7.57 x 10~ l6 cm 3 /s 
2=1.20 

Vi = 150 keV 
V e = 0.01 keV 
M— 1.68 

Injector efficiency = 0.80 
DC efficiency = 0.76 
Thermal efficiency = 0.50 

Unpolarized injection 
U ae = 2242 keV 
U ai = 416 keV 

Polarized injection 
U ae = 2681 keV 
U ai = 487 keV 

Alpha trapping enhancement = 1 .09 
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(x 0.01) 




T (keV) 

Fig. 11. Q contours for D-T reactor. 


(x 100) 


For a simple mirror fueled with polarized nuclei, 
both the <50% reaction index and the higher mean ion 
energy should enhance expected Q values. The ratio 
Qpoiarized / Qdassicai shown in Fig. 10 for various frac- 
tional burnup rates. For the LMRD parameters, a Q 
enhancement factor of 1 .63 is found here. The impact 
of this Q enhancement on nr requirements is shown in 
Figs. 11 and 12. 


DISCUSSION 
System Efficiency 

Here, Q is chosen as the reactor parameter of in- 
terest because it simply displays the breakeven Q = 1 
and because it can be written in terms of general phys- 
ics specifications. For operating power plants, however. 
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Fig. 12a. Effect of polarization on m requirements. 



(x 100) 


Ion Temperature (keV) 

Fig. 12b. Effect of alpha-particle heating on Q contours. 


the overall thermal efficiency may prove more useful. 
This overall efficiency is defined as the ratio of usable 
power to fusion power, or e = P net /P F - Decompos- 
ing the power terms as shown in the flow diagram of 
Fig. 13 and neglecting radiation and auxiliary power 
terms, the efficiency can be written in terms of Q as 


e 


Si 



+ 


( Pc/ P\ ) £ DC 



1 -Sit DC 

hQ 


where Pc/Pn is the ratio of total power in charged spe- 
cies to neutral species and & are the operating 

efficiencies of the thermal converter, the direct con- 
verter, and the injector, respectively. Substituting the 
LMRD mirror parameters, Fig. 14 summarizes the ef- 
fect of polarization on overall efficiency. 

Likewise, the recirculating power fraction, i.e., 
the ratio of power consumed to gross electrical pro- 
duced, is another important factor in mirror feasibil- 
ity evaluations and is intimately related to overall plant 
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Fig. 13. Reactor schematic. 



Fig. 14. Overall thermal efficiency versus Q. 
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Fig. 15. Fractional recirculating power versus Q. 


efficiency. Persistently high recirculating power require- 
ments alone account for 13% of the total capital out- 
lays estimated for the LMRD, and attempts to reduce 
this fraction are an important part of any cost reduc- 
tion scheme. This fraction can be written as 



and is shown in Fig. 15 for polarized and unpolarized 
injection and various Q values. 

Mirror Assessment 

Given that power cost estimates for the standard 
mirror reactor run four times greater than the Wisconsin 
tokamak estimates and ten times greater than current 
fission costs, the mirror’s future economic competitive- 
ness depends on Q enhancement schemes. As dem- 
onstrated in Fig. 16, modest improvements in Q can 
significantly reduce the cost of power, and a Q enhance- 
ment of 5 is expected to bring mirror costs into com- 
petition with other future energy sources such as solar, 
coal-fired magnetohydrodynamics, and the fast breeder 
in the range of 50 to 70 mil/kW -h (Ref. 19). 

While the Q enhancement of 1.63 found here for 
polarized nuclei falls short of these requirements, the 
potential for supplementing tandem or field-reversed 
mirror confinement with polarization holds some prom- 
ise. Estimates for these other confinement schemes 20 



Q/Qcl 

Fig. 16. Effect of Q enhancement on cost of power. 


are in the Q enhancement range of 3 to 5, and, coupled 
with a low-cost polarization program, may provide a 
mirror Q = 5 value or greater as required for future eco- 
nomic comparisons. 
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SUMMARY AND CONCLUSIONS 

The conceptual design of a mirror fusion reactor 
is analyzed to include polarized fuel. In discussing the 
various polarized designs, previous work has considered 
that the sin 2 1 ? angular distribution of fusion products 
might favor alpha-particle confinement and thus in- 
crease plasma heating. Compared with the classical un- 
polarized case, polarized fuel is found here to yield 
enhanced Q values (by a factor of 1.63). Estimated 
power is obtained from a Fokker-Planck representation 
that includes both direct and collisional losses. These 
improvements in alpha-particle confinement derive from 
polarized D-T nuclei that, in turn, produce a highly 
anisotropic alpha-particle source (sin 2 $) distribution. 
Following collisions, the alpha-particle products are 
tracked using a numerical solution and are found to 
maintain much of the high anisotropy initially fed to 
the reactor as polarized fuel. 

To extend the model into higher collective modes 
of plasma behavior, nonlinear instabilities associated 
with the loss cone should be considered. A host of 
anomalous sources for alpha-particle losses can be rep- 
resented, but nonlinear simulations may prove compu- 
tationally difficult without resorting to more powerful 
scaling of variables or spectral techniques beyond the 
scope of this paper. For polarized fuel, the principal 
outcome of enhanced heating may encourage further 
progress not only for simple mirror designs but also for 
more elaborate geometries such as tandem or sphero- 
mak reactors. 
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